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Abstract 

We pursue a geometrical approach to gravitational lcnsing theory. We present a 
survey of the background theory of General Relativity, including particular proper- 
ties of the Schwarzschild and Kerr solutions. Next we outline a proof of the Gauss 
Bonnet theorem and its applications to surfaces in optical geometry, as developed by 
G. W. Gibbons and M. C. Werner. Finally, we attempt to extend this geometrical 
approach to the axially symmetric Kerr spacetime, and arrive at an expression for 
the gravitational deflection angle in the equatorial plane. 



1 Introduction 



In 1916, Einstein published "The Foundation of the General Theory of Rel- 
ativity," an extension of his special theory of relativity to include gravity. From 
a contemporary perspective, physicists salute this achievement as the inaugura- 
tion of twentieth-century physics. Indeed, entire theories of scientific advancement 
as paradigm shifts may take the 1916 paper as their model: so Kuhn's notion of 
paradigm shift seems to stem from this impulse. Thus, the development of physics 
and the history of science derive from Einstein's radically new and intuitively sur- 
prising claims. In addition, the paper impelled new directions in experimentation. [5] 
Indeed, it was not long until a series of experiments verified his most radical thoughts. 
The theory came in a series of papers from 1911 to 1916, and was followed by the 
empirical research. 

Einstein's theory of gravity successfully reproduced Newtonian gravitation in 
domains with relatively weak gravitational fields. Moreover, it explained the extra 43 
arcseconds per century in the precession of Mercury's orbit, previously unaccounted 
for. One of the most striking predictions of Einstein's theory of General Relativity 
is the deflection of light due to the curvature of spacetime. In 1919, Eddington's 
observation and confirmation of this effect made Einstein an instant international 
celebrity. 

While Newton's theory of gravity established the idea of force, instead Ein- 
stein viewed mass itself as equivalent to the curvature of spacetime, in a precise 
sense. In contrast to all other forces, no known object is immune to the effects of 
gravity. Paths previously consider to be curved by Newton -say the orbit of the 
moon- were described by Einstein simply as straight lines in a curved spacetime, 
or geodesies. Therefore, while Newton held that light, composed as he believed of 
massless particles, would travel in a straight line through empty space, Einstein's 
theory implied the path of light, indeed of all objects, should curve around a mas- 
sive body. This is the gravitational lensing effect. In 1919 Sir Arthur Eddington, 
then secretary of the Royal Society and chief assistant to the Astronomer Royal of 
the Greenwhich observatory, received papers from Willem de Sitter describing Ein- 
stein's theory. He organized two expeditions to test the gravitational lensing effect 
by observing a solar eclipse. Despite Eddington's unfounded assumption that the 
data fall into one of three values, the expedition was considered the second piece of 
experimental support for Einstein's remarkable theory, following the calculation of 
the precession of Mercury in 1915. [3] 

The first known intergalactic lens, "Twin QSO," was discovered in 1979. [TT] 
Since then hundreds of microlensing events have been identified, spurring research 
and even allowing the detection of previously invisible celestial objects. Gravi- 
tational lensing theory studies geometrical optics in spacetime, connecting astro- 
physics, theoretical physics, and mathematics. Recent research develops a new 
method of analyzing lensing systems, one which connects the differential geome- 
try of these systems with their topology. [12] This paper will attempt to extend this 
work from spherically symmetric spacetimes to the Kerr spacetime. 
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2 The Schwarzchild spacetime 



In 1916, shortly after Einstein published his theory of G.R., Karl Schwarzschild 
found an exact solution to Einstein's equations that describes the unique structure of 
a nonrotating spherically symmetric spacetime. The Schwarzschild metric success- 
fully reproduces Newtonian gravity in the non-relativistic (weak field) limit. More- 
over, it allowed a calculation of the precession of the orbit of Mercury, accounting for 
the extra 43 arcseconds per century that had long been observed. The Schwarzschild 
metric also allowed the slight deflection of light by the sun to be calculated, and this 
was soon followed by Sir Arthur's expedition in 1919, which observed such deflection 
during a solar eclipse. This provided a strong confirmation of Einstein's theory of 
General Relativity, published just three years before. The components of the metric 
g are given in terms of the line element, ds 2 = g^dx^dx" 11 , as follows. 

ds * = _ ^ _ 2K) dt 2 + i (! _ *fy l dr 2 + gsin 2 (#)d0 2 + ^dO 2 

We choose units to set the speed of light c = 1, and we will follow the 
convention of using Latin indices to denote space coordinates, and Greek 
indices to denote spacetime coordinates. Since light is massless, light rays 
will follow null geodesies in accordance with G.R. Since the Schwarzschild 
metric is spherically symmetric, we may restrict ourselves to the equatorial 
plane without loss of generality, letting 6 = | and d9 = 0. Fermat's principle 
states that light rays follow a path of critical time. Furthermore, a spatial 
projection of the null geodesies will produce the light rays. [12] We find that 
the following metric describes the optical reference geometry: 

dt2 = V^7 dr2 + tAm)^ 2 

9mn — 9mn/(—9oo) 

Although our optical surface has only (2) spatial components, let us define 
in general a T^, 

v lv '■= \9 Xp (dx»g P v + d x»9 P v - d xP g^ u ) 

Where a sum is implied over p, in keeping with the Einstein summation 
convention, which we shall continue to use. This gamma essentially repre- 
sents the A - component of the covariant derivative of the derivative of the \x 
basis vector with respective to the v basis vector. To be precise: 

Next we calculate the Gauss Curvature K of this optical surface. 
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2.1 Calculation of the Gauss Curvature 

The Riemann curvature tensor R measures the failure of a vector to 
return to its original position following parallel transport about some in- 
finitesimal loop in a given plane. The tensor can be defined on any manifold 
Ai with an afline connection, and on a two dimensional manifold this ten- 
sor can be shown to have only one independent component. [7] Although the 
Gauss Curvature K can be defined as the product of the two principal cur- 
vatures k\ and &2, which are the eigenvalues of the second fundamental form 
of A4, Gauss's Theorema Egregium proves that K is in fact an intrinsic 
measure of curvature, and K can be expressed purely in terms of the met- 
ric. For a surface, it is not hard to show that the independent component 
of the Riemann tensor, say -R1212, is related to the Gauss Curvature K by 
K = R\2YilQi where g is the determinant of the metric. [7] We wrote the 
following Mathematica code to calculate the Gauss Curvature of the optical 
surface of the Schwarzschild spacetime. 



Clear [X, r, t, p, T, RAbcd, Rabcd] ; 
X ={r,0}; 

g = {{r A 2/(r - 2M) A 2, 0}, {0, r A 3/(r - 2M)}}; 
gi = Inverse[g]; 

r[A_,/i_,i/_]:=Sum [l/2gi[[A,p]] (d x[[u]] g[[p, n\\ + d x[M] g[[p, u]] - d x[[p]] g[[n, u]]) ,{p,2}] ; 
RAbcd[a_,b_, c_, d_]:=d x[[c]] T[a, d, b] - d x[[dl] T[a, c, b] + 

Sum[r[a, c, A]T[A, d, b], {A, 2}]-Sum[r[a, d, A]T[A, c, b], {A, 2}]; 
Rabcd[a_, b_, c_, d_]:=Sum[g[[a, z]]RAbcd[z, b, c, d], {z, 2}]; 
K[a., b_, c_, d_]:=Rabcd[a, b, c, d]/(g[[a, c]]g[[d, b]} - g[[a, d]}g[[c, &]]); 
FullSimplify[K[l,2,l,2]] 
We find that the Gauss Curvature is: 

_ M(3M-2r) 
r 4 

This expression varies only with the radius r, and the mass-like constant M, 
as we expect. 



4 



2.2 The geodesic equations 

Geodesic curvature is an intrinsic measure of the amount a curve is 
curving within a manifold. In General Relativity, all objects, in the absence 
of forces, follow a path of zero geodesic curvature (geodesies). Define a curve 
C : K — > M, such that p h- > x (p). The geodesic is then defined to be the 
non-normal component of the change in the curve C"s direction and speed: 

K g := < V x >x',x"/\x"\ >=g(V x ,x',x"/\x"\) 

From this definition and the one for T u \^, one arrives at the following differ- 
ential equation, which defines a geodesic curve x, parameterized by the affine 
parameter p\ . 

dx^ I x^x^_ Q 

dpf v\ dpi dpi 

We define the metric coefficients E, F, and G as follows: 

E:= < dr, dr >= 

F:= < dr, dep >= 

G:= < d<j), d<j) >= ^J^y 

Clearly, Ea, = = 0. This means our patch is Clairaut in r and so 
the geodesic equations for light, in the optical geometric surface of the 
Schwarzschild equatorial plane, reduce to: 

r" + §§ (r') 2 -ff (</>') 2 = 

0" + = o 

These differential equations do have an exact solution; however, there 
is a different method of deriving the asymptotic deflection and of light pass- 
ing by the Schwarzschild black hole, without having to solve or integrate 
the geodesic equations of motion. This method is due to [12], and man- 
ages to avoid having to deal with the curvature singularity of the black hole 
altogether. 
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3 Geometrical approach to lensing 



Now we will apply the method developed in [12] to find the asymp- 
totic deflection angle of a light ray, with impact parameter b, that passes by a 
Schwarzschild black hole of mass m = . We will utilize the Gauss-Bonnet 
theorem, which elegantly connects a surface's geometry to its topology. The 
Gauss-Bonnet theorem then states that given a 2-dimensional manifold Ai, 
with piecewise smooth boundary dAi, interior angles ij, geodesic curvature 
Kg, Gauss curvature K, and Euler Chi function x{M): 

/ / KdA + / K g ds + V (tt - ij) = 2n X (M) 

J JM JdM j 

Following [12] . we define the gravitational lensing system in the following 
way, and seek to find the asymptotic deflection angle of a light beam: 




Figure 1. 

Let S be the source of the light, O be the position of the observer, and 
L be the spherically symmetric mass. Let 71 be a light geodesic from S to 
O, and let b be the impact parameter of 71 with respect to L. Let 7 p be part 
of a sector of a circle centered at L and intersecting 71 at S and O at right 
angles. Let 4> he the angular coordinate and let r be the radial coordinate, 
both centered on L. Let 5 be the asymptotic deflection angle of 71. Let r$ 
and ro be the radial coordinates of S and O, respectively. We shall assume 
that b, rs, and ro m, so that the geodesic 71 is always in the weak field 
limit. Furthermore, we assume rs, ro 3> b, so that 71 may be approximated 
as the union of two straight segments. [12] 

Now let us find 5. Let D be the region bounded by 71 and 7 P . Let the 
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boundary dD be 71 traversed from S to O and 7 P from O to S. The angular 
range of 71 is tt + 5 and x(-D) = 1. We apply the Gauss-Bonnet theorem to 
find: 

/ f D KdA + k 9 ( 7p ) (l) d<\> + 2f = 2vr 

JT*«p(7p) ($)^-- = -/L™- 

We are interested in the asymptotic deflection angle, and so since the 
Schwarzschild metric is asymptotically flat, over 7 p we can take k 9 ( 7p ) % = 
1. Similarly, we can describe 71 as r (0) = b/sin{4>). [US p. 3] The equation 
then becomes: 

S = -Jf D KdA 

dA = y/det (g°pt)drd(j) = J ' {r ^ M yi drd(f) 
r* roo 2M { ^> drdcb 

JO Jb/sin{(t>) r ( 1 _m\ u '' a V 

~ 4M 
~ b ■ 

This matches the known formula for the deflection angle. See [12] for a 
fuller exposition of this geometrical approach to gravitational lensing. We 
will next examine the Gauss-Bonnet theorem in more detail, providing the 
main steps of a proof. 
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4 The Gauss-Bonnet theorem 



Let M be a 2-dimensional manifold with boundary dAi, and with 
interior angles ij. Then the Gauss-Bonnet theorem relates the integral of the 
Gauss curvature, plus the total angle curved out by the boundary of M to 
the Euler Chi function in the following way: 



K g ds+ / KdA + ^2(n - ij) =2ttx(M) 



dM J J M 



To prove this, first let a : [0, 1] — > dM be a simple closed unit-speed curve, 
parameterized by < s < 2tt. Let M. be parameterized by x(r,(f>), with 
X(u,v) G T X M. Let 

X u := d u X 
X v := d v X 

We can assume the parameterization has been chosen so that g (X u , X v ) = 0. 
Let 



E:=g(X u ,X u ) 
G:=g(X v ,X v ) 



Ei :- V 



y/\g(X»,X tt )\ 



E 2 := V 



y/\g(X v ,X v ) 



By construction, we see that g (E\,E2) = 0, and g (E\,Ei) = g (E2, E2) = 1. 
So we have an orthonormal basis. Define: 

w ij = ^a'EjEi 

W21 = V a iE\E2 

cos (6) :=g(a',E 1 ) 

a' = cos (<9) E 1 + sin (0) E 2 

V a ,a' = ((de)/(dt) + wax) (- sin (0) ^ + cos (0) E 2 ) 

K g := g (V a ,a', a') = (f + W21) (sin (6) 2 + cos (0) 2 ) = f + u; 21 
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oi = X u v! + X v v' 



Vo/tfi #2 = [(Si) u d«/ds + (^l)t, dv / ds \ E ? 

= ( [X UU /VE - X u EJ (2£ 3 / 2 )] f g + [X UV /VE - X u EJ (2£ 3 / 2 )] f ) ^ 
Note: 



G 11 



Ey 

2 



Therefore 

W2i = (2>/EG) £ + G tt / (2^^G) 



(in 



Recall the formula for the Gauss curvature K in terms of the metric: 
K = 

Now, 



/EG {{veg) v + (vSb) J 



d0 

W21 = S - S 



d0 



% = £„/ (2v^G) ^ - GJ (2VEG) 
Jf s ds-J Kg ds = f (E v / (2VEG) f s - GJ (2VEG) *) d S 
= J a EJ (2^EG) du - GJ (2^EG) dv 



Now, applying Green's theorem, we get: 



I Im 2 



Veg)J \ 



EG J v 



dudv. 



J j M K\J~EGdudv 
II M KdA 



So, 



/«M^ ds + /XM^ tW = /S * = *(2tt) -0(0) 



Furthermore, if we allow a piecewise-continuous boundary to our man- 
ifold, then we have corners which contribute to the total angle subtended, 
giving: 

f aM K g ds + / f M KdA + £\ [ir-ij) = 9 (2tt) - 0(0) = 2vr an integer. 



We shall argue that the integer multiplying 2ir above, call it xi-M), is 
indeed the Euler Chi function. Suppose we have some triangulation of Ai. 
Let T be the set of all triangles. We can sum the equation above for every 
triangle t G T. Note that each triangle is a topological disk. We therefore 
find that the sum of these equations is equal to 2irF, where F is the number 
of triangles. 



Et f d t K 9 ds + Jt KdA + ^j( 7T -itj) 



2ttF 



Note that the sum of the integrals of K over the triangles' areas is equal 
to the integral of K over the whole area of the manifold. Now note every 
interior edge of the triangulation is traversed once by a triangle in one di- 
rection, and once by the bordering triangle in the other direction. Therefore 
the two integrals of n g over such an interior edge cancel each other, and so 
their sum is zero; furthermore, if we remove such an interior edge, we have 
reduced the number of faces by one. Therefore we can change the sum of 
the integral of n g over the boundaries of the triangles, to an integral of n g 
over the boundary of the manifold, provided we also subtract 2tt times the 
number of interior edges from the right side of the equation. Now note that 
every interior vertex contributes a total of 2tt to the sum of the interior angles 
of the triangles. Therefore let us remove these interior vertices from the sum 
of the interior angles, and simply add 2tt times the number of interior ver- 
tices to the right hand side of the equation. To complete the proof, we note 
that the number of exterior edges must be equal to the number of exterior 
vertices, since the boundary of the manifold is a union of topological circles. 
Hence x(.M) = V — E + F, as desired. Therefore, we have found: 

Ibm K 9 ds +I Im KdA +T^ (t - ij) = 2n X (M). Q.E.D. 
Proof based from [3 p. 281-291]. 
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5 The Kerr Metric 



The Kerr metric is an exact solution of Einstein's equations which de- 
scribes the axially symmetric spacetime of a rotating spherically symmetric 
mass. It is a more realistic model of the spacetime of a black hole, as real 
black hole generally have non-zero angular momentum. 

An accurate formula for the asymptotic deflection angle of light around 
the Kerr black hole could provide a second direct method of measuring the 
angular momentum of black holes, and could thereby be used to test the idea 
of cosmic censorship, proposed by physicist Dr. Roger Penrose. |13j 

Cosmic censorship is the idea that the universe does not allow "naked" 
singularities. Singularities are the Achilles heal of Einstein's theory of Gen- 
eral Relativity since the theory breaks down there-the Riemann curvature 
tensor diverges. Physicists Dr. Penrose and Dr. Hawking have proposed 
that such singularities must be enclosed in a (black) event horizon to keep 
them from causally interacting with the rest of the universe. If this were not 
the case, eventually one would no longer be able to say what will happen 
anywhere in the universe after a certain time. |12j 

The Kerr metric is as follows: 

dg2 = _ A-aW(fl ) (fe2 _ Warsm*(e) dtd(f) + 

(a2+r2)2 -/ Asin2(g) sin 2 (9) d<t? + £dr* + p 2 d9 2 p. 309] 
A := a 2 - 2Mr + r 2 

p2 ._ r 2 _|_ fl 2 cog 2 ^ 

For the sake of simplicity, we shall restrict ourselves to the equatorial plane, 
so that 9 = w and d9 = 0. Because the metric is reflection-symmetric about 
9 = f, a trajectory which begins with p e = will continue to have p e = 0.[9l 
p. 359] Therefore a geodesic that begins in the equatorial plane will remain 
in it. The metric g becomes: 

ds 2 := g^dxW 

We again wish to find null geodesies, so we set ds 2 = 0, and solve for dt in 
order to find an optical metric. 

i^Ll dt 2 + iaM dtd(j)= ^ dr 2 

(V^2M-rdt + ^fc^) 2 = £<ir 2 +(r 4 + a 2 r 2 + 2a 2 Mr + ggg) 
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dt — y/ A(r-2M)^ r2 + JT^Mpdfi 2 (r-2M) 

This asymmetrical form of metric was first discovered and studied by Flanders. [8j 
We have found our optical Finsler metric for the optical reference geometry 
of the Kerr. Further, [2} gives several formulas for the Finsler equivalent of 
Gauss curvature, as well as a computational algorithm for computing it. 

5.1 Calculation of Gauss Curvature 

We hypothesize that the leading-order term in a = J 2 /M, of the 
optical deflection angle of the Kerr may come not from the asymmetrical 
Finsler part of the metric, but rather from the pseudo-Riemannian part. 
The following Mathematica code was written to arrive at an estimate for 
the Kerr deflection angle based on the pseudo-Riemannian (diagonal) part of 
the metric, ignoring the Finsler contribution, in order to test this hypothesis. 

Clear [X, r, t, fj,, T, RAbcd, Rabcd] ; 
X = {r,0}; 

A = r 2 - 2Mr + a 2 ; 

g = {{r 3 /(A(r - 2M)), 0}, {0, (r 2 a 2 + r 3 (r - 2M))/(r - 2M) 2 }}; 
gi = Inverse [g]; 

r[A_,/i_,i/_]:=Sum [l/2gi[[A,p]] (d x[[u]] g[[p, v}} + d x \\^g[\p, v]] ~ <9x[[ P ]]#[[/^ "]]) . {p, 2 I] I 

RAbcd[a_,b_, c_, d_]:=d x i [c]] T[a, d, b] - d X [[ d ]]^[a, c, b] + 

Sum[r[a, c, A]T[A, d, b], {A, 2}]-Sum[r[a, d, A]T[A, c, b], {A, 2}]; 

Rabcd[a_, b_, c_, d_]:=Sum[p[[a, z]]RAbcd[z, 6, c, d], {z, 2}]; 

K[a., b_, c_, d_]:=Rabcd[a, b, c, d]/(g[[a, c]]g[[d, b]} - g[[a, d]]g[[c, &]]); 

FullSimplify[K[l,2,l,2]] 

We arrive at the following expression for the Gauss Curvature K of our re- 
duced metric: 
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M(6a 2 (-M+r)+r(6M 2 -7Mr+2r 2 )) 
"~ = (2M-r)r> 

We find that the limit as a — >■ of this K reduces to the value found for the 
Gauss K of the optical geometry of Schwarzschild, as expected. 

5.2 Extension of the Geometric Approach 

Now we attempt to extend the geometric approach which we applied 
to the Schwarzschild spacetime to the Kerr spacetime. Starting from Gauss- 
Bonnet and using the expression for K above, we find: 

dA = y^\drdcP = ^( a^m) ) ( (r "2M) 2 ) drc ^ 
6 = ~ I Id K (^§^372 drd<f> 

~ ~ Jo Jb/Sin[ct>] K (r-2M)3/? drd( t ) 

_ r-TT poo Mr 5 / 2 (6a 2 (r-M)+2r 3 (l-2M/r)(l-3M/2r)) , , , 

— Jo Jb/Sia[<j>] r 15 /2(l-2M/r) 5 /2 ara f 

= Jo /vSin M M 6a 2 (l - M/r)/{r\l - 2M /r)^)drd<t> + 

Jo JyanM 2M (! " 2M / r )(! " 3M/2r)/(r 2 (l - 2M/r) 5 / 2 )rfr# 

- Jo JvSinM ^ 2 Mr-*drd<f> + ft J^ m 2Mr^drd^ 
= Jo / 6 ?Sin W 6a 2 Mr-4dr# + AM/b 

= 2a 2 M/b 3 £ Sin[</>] 3 d(f> + 4M/b 
= 2Ma 2 /6 3 (4/3) +AM/b 

= ¥(i + !S) ° 

And so we have found the contribution to the deflection angle from the 
pseudo-Riemannian part of the metric shows up at second order in a. 
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6 Discussion 



In [I] and [TO], expressions are derived for the deflection angle of light 
curving around a Kerr blackhole, in the weak-deflection limit, to third order 
in a and in M, where M is the scaled mass of the blackhole, and a = jj 
is a measure of the angular momentum. Let b be the impact parameter. 
Restricting their expression to the equatorial plane, the angle of deflection 
for a retrograde photon reduces to: 

5= « + (f 7r M 2 + 4Ma)/b 2 + {±fM 3 + WirM 2 a + 2Ma 2 )/b 3 .[10, p. 5] 

We note the second order term in a, 2j£ a 2 , is reasonably close to the 
answer we got from the Riemannian part of the metric: | jjf a 2 ; however, the 
first order terms in a do not come from the Riemannian part of the metric, 
but rather from the Finsler geometry. Indeed, while the diagonal part of the 
metric does not distinguish between a positive displacement in <fi and a nega- 
tive one, in reality we should find that retrograde light rays, moving against 
the rotation of the blackhole, should have an increased deflection angle, while 
prograde light rays should have a decreased deflection angle. This is caused 
by the light ray spending more time near the blackhole in the retrograde case, 
and less time close to the black hole in the prograde case. [10] If we define 
to equal +1 for a retrograde ray, and —1 for a prograde one, the expression 
for the magnitude of the asymptotic deflection angle from [10J becomes: 



AM 



+ 



15tt M 2 



+ 



128 M 3 
3 ~W 



+ 5j 



AM 



+ 



IOttM 2 



So we see that retrograde light rays have a larger deflection angle than 
prograde ones, and that while the terms linear in a depend on which side of 
the blackhole the light ray passes, the a 2 term has no such dependence. 
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7 Conclusion 



The curvature of the Kerr and Schwarzschild spacetimes is negative. 
This should imply that locally parallel geodesies diverge. The topology of 
these spacetimes is therefore essential to the convergence of light rays, and 
hence to the gravitational lensing effect. The Gauss-Bonnet theorem, which 
relates the local geometry of a two-dimensional manifold to its global topolog- 
ical structure, directly connects the global properties of the optical geometry, 
necessary for gravitational lensing, to the measurable deflection angles. The 
theorem thereby provides an effective, informative framework for examining 
lensing. 

In conclusion, we have found that the Gauss curvature derived from 
the pseudo-Riemannian part of the Kerr metric does not account for the 
first-order effect of the angular momentum quantity a = jj , but instead ap- 
pears at second order in a. Perhaps future work could examine whether the 
Gaussian curvature analogue in Finsler geometry can be used in a manner 
similar to the Gauss curvature, to deduce a measure of the asymptotic optical 
deflection angle in the Kerr spacetime. 
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